Self-force on extreme mass ratio inspirals via curved spacetime effective field theory 
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In this series we construct an effective field theory (EFT) in curved spacetime to study gravi- 
tational radiation and backreaction effects. We begin in this paper with a derivation of the self- 
force on a compact object moving in the background spacetime of a supermassive black hole. The 
EFT approach utilizes the disparity between two length scales, which in this problem are the size 
of the compact object r m and the radius of curvature of the background spacetime 1Z such that 
e = r m /TZ -C 1, to treat the orbital dynamics of the compact object, described as an effective point 
particle, separately from its tidal deformations. The equation of motion of an effective relativistic 
point particle coupled to the gravitational waves generated by its motion in a curved background 
spacetime can be derived without making a slow motion or weak field approximation, as was assumed 
in earlier EFT treatment of post-Newtonian binaries. Ultraviolet divergences are regularized using 
Hadamard's partie finie to isolate the non-local finite part from the quasi-local divergent part. The 
latter is constructed from a momentum space representation for the graviton retarded propagator 
and is evaluated using dimensional regularization in which only logarithmic divergences are relevant 
for renormalizing the parameters of the theory. As a first important application of this framework 
we explicitly derive the first order self-force given by Mino, Sasaki, Tanaka, Quinn and Wald. Going 
beyond the point particle approximation, to account for the finite size of the object, we demonstrate 
that for extreme mass ratio inspirals the motion of a compact object is affected by tidally induced 
moments at 0(e 4 ), in the form of an Effacement Principle. The relatively large radius-to-mass ratio 
of a white dwarf star allows for these effects to be enhanced until the white dwarf becomes tidally 
disrupted, a potentially 0(e 2 ) process, or plunges into the supermassive black hole. This work pro- 
vides a new foundation for further exploration of higher order self force corrections, gravitational 
radiation and spinning compact objects. 



I. INTRODUCTION AND MAIN POINTS 

In two previous papers [TJ O [3], using a stochastic 
field theory approach based on open system concepts, 
we derive the scalar, electromagnetic and gravitational 
self-force to leading order on a particle moving in an ar- 
bitrary curved background spacetime. We begin with 
the particle following a quantum mechanical path while 
interacting with a linear quantum field (HOE]. The 
conditions on a stochastic field theory (for open systems) 
to emerge from a quantum field theory (of closed sys- 
tems) are that the mass and size of the particle are large 
enough so the particle worldline is sufficiently decohered 
from its interactions with the quantum fluctuations of 
the field that it can be considered as quasi-classical, and 
yet sufficiently small that quantum fluctuations manifest 
as classical stochastic forces [7] . 



A. Quantum, Stochastic and Effective Field 
Theories 

When there is a significant discrepancy between the 
two mass (or energy or length) scales in a problem, as 
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in the extreme mass ratio binary systems under con- 
sideration here, one could use an open system stochas- 
tic description for their dynamics, such as developed in 
[H [2] • When this discrepancy is huge (such as between 
the QCD and GUT scales in particle physics) the stochas- 
tic component is strongly suppressed in the wide range 
between the two scales, away from the threshold region 
|8J . Then the stochastic field theory description will give 
rise to an effective field theory (EFT) description [3] for 
the motion of the small mass subsystem. Due to the 
large separation in the mass scales quantum loop correc- 
tions from the field and the intrinsic quantum mechan- 
ical worldline fluctuations are very strongly suppressed. 
These two factors render a quantum field theory (QFT) 
into a stochastic field theory (with sufficiently decohered 
histories) and in turn (with sufficiently small stochastic- 
ity) an effective field theory for the dynamics of the re- 
duced systems. We shall explain the essence and demon- 
strate the advantages in taking a field theory approach 
to treat radiation-reaction of classical systems. The ap- 
plication of EFT to the treatment of gravitational ra- 
diation from post-Newtonian (PN) binary systems was 
first introduced by Goldberger and Rothstein [ID] . Our 
formulation of a curved spacetime effective field theory 
(CS-EFT) goes beyond with two important features: it 
is for any curved spacetime background and there is no 
slow motion or weak field restrictions [69] . 

In this and three subsequent papers [IT] [12] [13] we 
construct a CS-EFT and apply it to derive the self-force 
on a compact object moving in an arbitrary curved back- 
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ground For concreteness, the background spacetime is as- 
sumed to be that of a supermassive black hole (SMBH) 
with curvature scale TZ much larger than the size of the 
compact object r m . The smallness of the ratio e = r m /lZ 
makes it a good expansion parameter for a perturbation 
theory treatment describing the extreme mass ratio inspi- 
ral (EMRI) of the compact object. These binary systems 
are expected to be good candidates for detecting grav- 
itational wave signatures using the space-based gravita- 
tional wave interferometer LISA [14] . We now give some 
realistic numbers for astrophysical processes in this cat- 
egory to delimit the range of validity for carrying out 
these perturbation expansions for these EMRI sources 
detectable in LISA's bandwidth . However, we emphasize 
that the formalism we construct here is of a sufficiently 
general nature that it can be applied to any compact ob- 
ject moving in an arbitrary curved background, including 
those spacetimes sourced by some form of stress-energy 
and those possessing a cosmological constant. 



B. Relevant scales in EMRIs 

Consider the motion of a compact object (a black hole, 
neutron star or white dwarf with a mass m ranging from 
about 1 to 100 solar masses) moving through the space- 
time of a SMBH with a mass M ~ 1O 5 ~ 7 M . We have in 
mind that the compact object moves in a stationary back- 
ground provided by the supermassive black hole, such as 
the Schwarzschild or Kerr spacetimes. Such spacetimes 
are appropriate for a description of the EMRI in which 
the compact object is bound by the gravitational pull of 
the SMBH. By emitting gravitational waves the binary 
system loses energy until the compact object plunges into 
the SMBH. The emission of gravitational radiation from 
such a system is expected to be detected with the antic- 
ipated construction and launch of the LISA space-based 
interferometer [14] . 

It is believed that most SMBHs lurking in the middle 
of galaxies, which are thought to host the prime sources 
of gravitational wave emissions detectable by LISA, are 
spinning and clean in the sense that most, if not all, of 
the surrounding material has already fallen into the black 
hole. (Active galactic nuclei are a notable exception [T5].) 
Because of this the Kerr background is perhaps the most 
astrophysically relevant spacetime for the extreme mass 
ratio inspiral. The Kerr solution is vacuous {R^ v = 0), 
stationary and stable under small perturbations [TB] and 
possesses two Killing fields. The first Killing field (£ Q ) 
is time-like and describes time-translation invariance ev- 
erywhere outside of the ergoregion. The second (ip a ) is 
space-like and describes the axial symmetry of the space- 
time. 

The Ernst [T7] and Preston-Poisson [TH] geometries de- 
scribe a black hole immersed in an external magnetic 
field. From an astrophysical viewpoint, the external mag- 
netic fields that a black hole at the center of a galaxy ex- 
periences are relatively weak and unlikely to significantly 



affect the motion of the compact object until a very high 
order in the perturbation theory. 

There are two relevant length scales in EMRIs. The 
smaller scale is set by the size of the compact object 
itself, denoted r m . For an astrophysical black hole its 
radius is r^h = ^G^ra ~ m/m^ where — 32ttGn 
in units where h = c = 1 |70j . For a neutron star with 
a mass « IAM@ and a radius of 10 — 16 km it follows 
that r ns « 4.8 — 7.7 'Gntti ~ m/rnHj. Therefore, it is to 
be expected that the size of the compact object, be it a 
black hole or a neutron star, is of the order of its mass 

cza - 

The second relevant scale is the radius of curvature 
of the background spacetime, TZ. We take TZ to be the 
following curvature invariant 

n=(R^R^y 1/A . (i.i) 

For a (possibly rotating) stationary SMBH the radius of 
curvature is 

/ rni,r 3 

where r is the typical orbital distance for the compact 
object away from the central black hole. For example, r 
is the geometric mean of the semi-major and semi-minor 
axes of a compact object in an inclined elliptical orbit. 
In an approximately circular orbit r is the orbital radius 
and for a particle moving faster than the escape velocity 
r is the impact parameter. 

In the strong field regime where r ~ M/m^ the cur- 
vature scale is also ~ M/m^ implying that r/lZ ~ m/M 
whence a perturbative expansion in e is equivalent to one 
in m/M. 

The typical variation in time and space of the back- 
ground is <^ TZ. The wavelength A of radiated metric 
perturbations from the compact object in a bound orbit 
is 

/m 2 ,r 3 

which shows that the wavelength of the gravitational 
waves does not provide a separate scale independently 
from TZ. 



C. The CS-EFT approach: Issues and main 
features 

The effective field theory approach was first introduced 
for the consideration of gravitational radiation from post- 
Newtonian binary systems in |10j . spinning compact ob- 
jects in [THl HOI [2D 122, and dissipative effects due to 
the absorption of gravitational waves in [SJJ HI]- See 
[25j[26l[27] for introductory reviews. Let us call these the- 
ories PN-EFT: they are constructed to describe the mo- 
tion of two slowly moving compact objects in a flat back- 
ground. In particular, the compact objects are treated 
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as effective point particles, the worldlines of which carry 
non-minimal operators describing the multipole moments 
from companion-induced tidal deformations as well as 
possible spin degrees of freedom and other intrinsic mo- 
ments. Below we describe briefly some general features 
of EFT and the specific differences between our new CS- 
EFT approach and the existing PN-EFT. 

The use of point particles to source the metric per- 
turbations about the flat background spacetime (note 
that the high frequency waves in a quantum descrip- 
tion corresponds to massless spin-two particles, the gravi- 
tons, in flat space quantum field theory) prompts the 
appearance of divergences. Fortunately there exists a 
well-established bank of tools and techniques in quan- 
tum field theory for regularizing these divergences and 
rcnormalizing the parameters and coupling constants of 
the theory. A theory is renormalizable in the effective 
field theory sense if observables are calculated in the low 
energy limit: the divergences can always be absorbed 
into a renormalization of the coupling constants of the 
infinite number of non-minimal worldline operators. The 
use of dimensional rcgularization is particularly useful in 
effective field theories because the renormalization group 
equations are mass-independent for this scheme indicat- 
ing that only logarithmic divergences contribute to the 
renormalization procedure |72], 

Our CS-EFT approach differs from this group of work 
in two ways. First, we work with an arbitrary curved 
spacetime. Second, we allow for the compact object to 
move with relativistic speeds in strong field regions of 
the background spacetime. The post-Newtonian effective 
field theory of fTUj treats bodies moving slowly through 
a weak gravitational field. 



1. In-In formulation for real and causal equations of 
motion 



Technically there are also fundamental differences. To 
derive real and causal equations of motion we empha- 
size the need to use the closed-time-path (CTP) inte- 
gral formalism based on an in-in generating functional 
[SSI 111 |13 [311 [31 133] ('in' and 'out' here refer to the ini- 
tial and final vacua used to define the vacuum transition 
amplitude in the generating functional). The authors of 
[TU] use the in-out formalism which is acceptable for field 
theories in a flat background spacetime since the in- and 
out- vacua are equivalent up to an irrelevant phase in that 
special case. In the presence of spacetime curvature, the 
difference becomes serious, as we will show in for 
the EMRI scenario. The in-out formalism can calculate 
matrix elements in scattering processes, but not expecta- 
tion values of physical observables in real-time evolution. 
The equations of motion for the effective particle dynam- 
ics obtained from an in-out formulation are not generally 
causal. An initial value formulation of quantum field the- 
ory via the CTP generating functional is the only correct 
way for the description of the system's evolution - it guar- 



antees real and causal equations of motion for the particle 
dynamics [3T| 133]. 



2. Effective point particle description 

Another important ingredient in our construction is 
an effective point particle description for the motion of 
the compact object. Going beyond the point particle ap- 
proximation is necessary to include the effects of tidal 
deformations induced by the background curvature as 
well as the effects from spin and other intrinsic moments. 
Following [10 , we introduce all possible terms into the 
point particle action that are consistent with general 
coordinate invariance and reparameterization invariance 
(and invariance under SO(3) rotations for a non-spinning 
spherically symmetric compact object). By implement- 
ing a matching procedure using coordinate invariant ob- 
servables we can match the observables of the effective 
point particle theory with the long wavelength limit of 
observables in the full "microscopic" theory to determine 
the values of the coupling constants of the non-minimal 
terms. As we will show in Section II VI this allows us to 
deduce the order at which finite size effects affect the 
motion of the compact object through the statement of 
an Effacement Principle. To our knowledge this has not 
been explicitly given in the literature before for the EMRI 
scenario. 



3. The power counting rules 

Power counting is a generalization of dimensional anal- 
ysis. In our perturbative treatment it is crucial for deter- 
mining how the Feynman rules scale with the parameter 
e. Once the scaling of the Feynman rules are known we 
determine all of the tree-level Feynman diagrams that 
appear at a particular order. Those diagrams contain- 
ing graviton loops are safely ignored. We also assemble 
the diagrams that include the non-minimal worldline op- 
erators describing the finite size of the compact object. 
Significantly, this allows us to determine the order in e 
at which finite size effects enter the particle equations 
of motion. With the power counting rules the EFT ap- 
proach becomes an efficient and systematic framework 
for calculating the self-force to any order in perturbation 
theory. Furthermore, by knowing how each Feynman di- 
agram scales with e we can study a particular physical 
interaction that is of interest by focusing our attention on 
a single diagram or on a few diagrams without having to 
calculate every contribution that appears at that order 
and at lower orders. For example, the leading order spin- 
spin interaction (spin here refers to that of the compact 
object, not of the SMBH) contributes to the self-force at 
third order in e for a maximally rotating compact object 
and can be calculated from the appropriate Feynman di- 
agram [3] [T3]. The power counting rules, the Feynman 
rules and their scaling with s are derived in Sections III B 
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4- Divergences and Regularization 

In Section |IIID| we propose a method for regularizing 
the divergences that appear in the effective action. Our 
approach utilizes a mixture of distributional and momen- 
tum space techniques within the context of dimensional 
regularization. We know from previous work that the 
finite part of the self- force is generally non-local and his- 
tory dependent. However, the ultraviolet divergences are 
quasi-local and independent of the history of the effective 
point particle's motion. To isolate the quasi-local diver- 
gence from the non-local finite part we use the method of 
Hadamard's partie finie, or finite part, from distribution 
theory. (See Appendix [B] for a brief review of the defini- 
tions and concepts of distribution theory relevant in this 
work.) 

After isolating the divergence from the non-local, finite 
remainder we then use a momentum space representation 
for the propagator in a curved background, first derived 
for a scalar field by Bunch and Parker [34], to calculate 
the divergent contributions. Their method is straight- 
forward but not efficient for higher spin fields, including 
gravitons in a curved space. (See also the related work 
of [35].) 

A novel method applicable for any tensor field is de- 
veloped in |36j for computing the momentum space rep- 
resentation of the Feynman propagator. The method 
is sufficiently general to do the same for any quantum 
two-point function, including the retarded propagator 
D™p^, s ,(x,x'). This approach makes use of diagram- 
matic techniques borrowed from perturbative quantum 
field theory. In Riemann normal coordinates, we expand 
the field action in terms of the displacement from the 
point x. The series can be represented in terms of Feyn- 
man diagrams, which allows for an efficient evaluation of 
each term in the expansion. Furthermore, we prove that 
some of the diagrams are zero to all orders. This iden- 
tity is not recognized in |34j even though its relation to 
certain steps made in their calculations is evident. 



D. MST-QW Equation 

Assembling all these essential ingredients, in Section 
|III| we give a demonstration of how the curved spacetime 
effective field theory construction is implemented, outline 
the steps in the regularization of divergences, perform an 
actual calculation of the effective action and from it de- 
rive the equation of motion for the compact object in- 
cluding the effect of gravitational self- force. As an appli- 
cation we work to first order in the mass ratio (i.e., e) and 
obtain the well-known Mino-Sasaki-Tanaka-Quinn-Wald 
(MST-QW) equation [37l[38]. This also sets the stage for 
calculating the second order self- force, the emitted grav- 
itational waves and the motion of compact objects with 



spin [njHU [Si- 
ll. EFFECTIVE FIELD THEORY FOR 
POST-NEWTONIAN BINARIES AND EMRIS 

Before proceeding to construct a curved spacetime ef- 
fective field theory for extreme mass ratio inspirals we 
briefly summarize the original work of |10j . which intro- 
duces effective field theory techniques for describing post- 
Newtonian binary sources of gravitational radiation. 

A. EFT in flat space for post-Newtonian binaries 

The aim of [10] is to describe the motion of two slowly 
moving bodies through a weak gravitational field using 
effective field theory techniques in order to generate a 
perturbative expansion in powers of the relative velocity. 
One of the many benefits of using an effective field theory 
approach is that the method is systematic and efficient 
so that there is in principle no obstacle calculating to any 
order in the velocity. 

The authors in [10] start by replacing the compact ob- 
jects with effective point particles. These are described 
by an action consisting of the usual point particle action 
plus all possible self-interaction terms that are consis- 
tent with general coordinate invariance and reparamc- 
tcrization invariance of the worldline. Then, the in-out 
generating functional is introduced to derive an effective 
action 

exp ^iS ef f[z] | 

= J Vh^ exp ^iS pp [z,ri + h] + iS[r] + h]^ (2.1) 

where S pp is the effective point particle action, S[ri + h] is 
the Einstein- Hilbcrt action for the full spacetime metric 
and z a are the coordinates of the particle worldline. 

Before integrating out the metric perturbations the au- 
thors observe that it is useful to separate into poten- 
tial H^ v and radiation h^ u contributions 

hfiu = + h^y. (2.2) 

This is suggested by the fact that the slowly moving bod- 
ies see a nearly instantaneous gravitational potential and 
yet radiate gravitational waves due to their relative accel- 
erations. This decomposition is also required to make the 
Feynman diagrams all scale homogeneously with the rel- 
ative velocity, v. In this way, the perturbative expansion 
in v is consistent and can be constructed to any order. 

Integrating out the potential gravitons using pertur- 
bation theory yields a theory of point particles moving 
in potentials. The radiation gravitons and the particle 
worldlincs are non-dynamical at this stage and can be 
treated as external sources. In this effective theory, valid 
at the orbital scale of the binary, the authors derive the 
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Einstein- Infcld-Hoffman potential as a check of their 
method. 

The last effective theory the authors have constructed 
involves integrating out the radiation gravitons. They 
then derive the famous power spectrum for quadrupo- 
lar gravitational radiation by calculating the first non- 
vanishing contribution to the imaginary part of the effec- 
tive action; the real part of the effective action generates 
equations of motion while the imaginary part is related 
to the power of the emitted gravitational radiation. 

Using an effective field theory approach it is not too 
surprising that some of the parameters of the theory un- 
dergo classical renormalization group (RG) scaling. In 
fact, the appearance of such RG scaling is used by the 
authors to show that there are no finite size effects up 
to v e order. In their words, "whenever one encounters 
a log divergent integral at order v 6 in the potential, one 
may simply set it to zero. Its value cannot affect physical 
predictions." [TU] This is how they resolve the problem of 
the undetermined regularization parameters that appear, 
at third post-Newtonian order, from regularizing the sin- 
gular integrals encountered with the traditional PN ex- 
pansion techniques. See [40 and references therein for a 
complete discussion of the regularization ambiguity. 



B. EFT in curved spacetime for EMRIs 

Our construction of an EFT does not rely on the slow 
motion of the bodies nor on the assumption that they 
move through a weakly curved region of spacetime. Quite 
the contrary, we allow for the compact object to move rel- 
ativistically through the strong field region of the SMBH 
background spacetime. As a result, the metric pertur- 
bations generated by the motion of the compact object 
cannot be partitioned simply into an instantaneous po- 
tential and radiation contributions. 

Within LISA's observable frequency bandwidth, e 
takes values between ~ 1CP 5 and 1CP 7 . Being so small 
over the dynamical time scale of the inspiral implies that 
e is a good parameter for building a perturbation theory 
within the context of effective field theory. 

Utilizing the dissimilar magnitudes of the compact ob- 
ject's size and the background curvature scale, we can 
construct two kinds of effective field theories. The first 
describes the compact object, in isolation from other ex- 
ternal sources, as an effective point particle. By allowing 
for all possible worldlinc self-interaction terms that are 
consistent with the symmetries of the theory we can ac- 
count for the tidal deformations, spin and intrinsic mul- 
tipole moments that the compact object may experience 
when it does interact with external sources. The second 
EFT is valid at scales ^> 1Z and results from integrat- 
ing out the metric perturbations. The resulting theory is 
that of an effective point particle subjected to a self-force 
from the gravitational radiation reaction as a result of its 
motion in the background spacetime, with the force and 
the radiation evolving self-consistently. Using a match- 



ing procedure we can establish the values of the coupling 
constants appearing in the effective point particle action. 

C. EFT of a compact object 

In applying the EFT formalism we first construct an 
effective point particle theory for the small mass m. This 
allows for a point particle description of the compact ob- 
ject's motion through the background spacetime while 
taking into account any tidally induced moments, or fi- 
nite size effects, that might affect its motion. An example 
is provided in classical electromagnetism wherein long- 
wavelength radiation scatters off a small, metallic sphere 
possessing no residual charge. The interaction of the field 
induces a time-varying dipole moment on the surface of 
the sphere. Since the wavelength is much larger than the 
sphere, this system can effectively be described by radia- 
tion interacting with a point particle carrying an induced 
dipole moment on its worldline. 

In the full theory describing the motion of a neutron 
star and the dynamics of the spacetime metric it moves 
in, the total action is given by 

3t„t = S[g] + S ns [g;p,p...]. (2.3) 

The quantities in the neutron star action, p,p, are 
the appropriate hydrodynamic variables necessary to de- 
scribe the internal dynamics of the neutron star whatever 
its equation of state. If the compact object under con- 
sideration is a small black hole then there is only the 
dynamics of vacuum spacetime to consider, which is de- 
scribed entirely by the Einstein-Hilbert action 

Slg} = 2m 2 pl Jd 4 xg 1 / 2 R (2.4) 

where R is the Ricci curvature scalar of the spacetime 
and g is the absolute value of the metric 's determinant. 

The effective point particle description of the com- 
pact object is constructed by "integrating out" the short 
wavelength gravitational perturbations up to the scale 
r m ~ m/rript. In doing so we introduce an effective point 
particle action S pp to describe the motion of the compact 
object moving in the combined geometry of the back- 
ground SMBH plus long wavelength gravitational per- 
turbations so that the total action becomes 

Stot = S[g] + Spp[z, g]. (2.5) 

Here z a (X) are the coordinates of the particle worldlinc 
and A is its affine parameterization. 

Being a description of the extended compact object the 
effective point particle action should include all possible 
self-interaction terms that are consistent with the sym- 
metries of the theory [5], which are general coordinate 
invariance and worldlinc rcparamctcrization invariance. 
For the discussion here, we will assume that the compact 
object is perfectly spherical when removed from external 
influences (e.g. background curvature) so that it carries 
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no permanent moments. This implies, for example, ex- 
cluding spinning compact objects in our construction, at 
least for now [73]. Hence, S pp should also be invariant 
under SO(3) rotations. From these considerations the 
most general such action is 



z,g] = —to J dr + crJ dr R+ cy J 

+c E fdrE^E^ + c B JdrB^B^ 



(2.6) 



which is effectively an expansion in powers of the compact 
body's radius r m over the much longer wavelength of the 
gravitational waves A. This can be interpreted as a mul- 
tipole expansion where the multipoles carry information 
about the induced moments that the background cur- 
vature and long wavelength gravitational perturbations 
impart to the compact object. We showed earlier that in 
the strong field region of a SMBH the wavelength of the 
metric perturbations is the same order as the curvature 
scale of the background spacetime, A ~ 1Z. This implies 
that the terms in the above multipole expansion have a 
definite scaling with e, which we will later demonstrate 
in Section ITVl 

The symmetric and traceless tensors £ M „ and Bp_ v are 
the electric and magnetic parts of the Weyl curvature, 
defined as 



f 



^ pavp <* 

a/3 ->^LP 
up 



£p.a@\C r '„„z : 



(2.7) 
(2.8) 



where z a is the particle's 4- velocity. Recall that these 
tensors are orthogonal to z a and C^ au /3 — Rp av p for a 
Ricci flat spacetime. 



The terms in the effective point particle action (2.6) 



proportional to the Ricci curvature vanish at leading or- 
der in e . The equations of motion for the full metric 

Rpu ~ \g^R = TPP(z,g) = 0(e), (2.9) 

where Tj^ is the stress-energy tensor of the effective point 
par ticle , can be used to set the second and third terms 
in (2.6 1 to zero at leading order. Equivalently, for the 



term proportional to cr, for example, we can redefine 
the metric g^ lL , in terms of a new metric g' through the 
field redefinition [10] 



SVO) = g'^ u (x) 



2m^ 



dr 



5 4 {x-z{t)) 



g 



11/2 



.(2.10) 



This transformation implies that the Einstein-Hilbert ac- 
tion is, to linear order in the (small) arbitrary parameter 

2m 2 pl J dtxgWRijg) 

= 2m 2 pl J d i xg /1 / 2 R(g , ) + ^ J dr R. (2.11) 



The term in S pp linear in R then appears with the con- 
stant cr + £, which can be set to zero since £ is arbitrary. 
One can find a similar field redefinition to remove the 
term proportional to cy. (See [H] H2] for a more thor- 
ough discussion of removing such terms from the effective 
field theory.) Using the metric field equations or, equiv- 
alently, performing a field redefinition of the metric one 
can remove all occurrences of the Ricci tensor in the effec- 
tive point particle action. It follows that the non-minimal 
couplings in S pp contain terms that depend only on the 
Ricmann curvature tensor. 

These field redefinitions allow for the effective point 
particle action (2.6 1 to be written as 



Spp[z,g] = -to J dr + c E J drE^E^ 

+c B f drB^B^ + ■■■ . (2.12) 



In a later section we will show that the non-minimal cou- 
plings in S pp are entirely negligible for calculating the 
MST-QW (first order) self-force. The MST-QW self- 
force is sufficient for computing gravitational waveforms 
and generating templates for LISA to detect gravitational 
waves from EMRIs [71]. However, a more precise deter- 
mination of the masses, spins, etc. of the binary con- 
stituents requires more accurate higher order templates, 
which can be computed by knowing the higher order con- 
tributions to the self- force [4"5l I44j . 

In the next section we derive the first order equations 
of motion for the compact object using the CS-EFT ap- 
proach. These equations, which describe the self-force on 
the mass to, were previously found by Mino, Sasaki and 
Tanaka |37j using matched asymptotic expansions and 
independently by Quinn and Wald [38] using axiomatic 
methods. In principle, we can compute the formal equa- 
tions of motion to higher orders in e thereby extending 
the work of [37J and [3H] . The second paper in this series 
[11 will give results through the second order in e . 



III. CS-EFT DERIVATION OF MST-QW 
EQUATION FOR FIRST ORDER SELF-FORCE 

In the previous section, we outlined the construction of 
an effective field theory that replaces the extended com- 
pact object by an effective point particle. In this section 
we construct an EFT for the motion of the effective par- 
ticle by integrating out the metric perturbations at the 
scale of the radius of curvature 1Z. In doing so, we derive 
the MST-QW self-force on the compact object. 

Denote the background (unperturbed) metric by g^ v so 
that the total metric is given by the background plus the 
perturbations generated by the presence of the moving 
compact object 



full 

7 fll/ 



rripi 



(3.1) 
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The metric perturbations /i M „ are presumed to be small 
so that \h^ u \ <C rripi. We will occasionally make use of 
the shorthand notation 



7 _ 'V 

tL l±V 



K 

rripi 



(3.2) 



for the (dimensionless) ratio of the metric perturbation to 
the Planck mass. From (2.5 1 the total action describing 



the interactions between the metric perturbations and 
the particle is given by the sum of the Einstein-Hilbert 
and effective point particle actions, 

S to t[g + h,z] = S[g + h] + S pp [g +h,z]. (3.3) 

We expand the Einstein-Hilbert action in orders of 

1 



S[g + h] = 



x g 



2h al 3, 1 h aT ' f} - h aP . a h a ^ 



2h- r 



V 



1 



rf 4 a;a 1/2 V; Q V" Q2 " plP2 



(n) 

n— 3 

x ^ p 1 h aia2 ^/ p 2 h a3a4 h a5aG • ' • h a2 

5 (2) + 5(3) + . 



-lQ2n 

(3.4) 



where denotes the part of the action containing 

terms proportional to n factors of /i M „. The quadratic 
contribution is the kinetic term for h^ v and provides the 
propagator corresponding to some appropriate boundary 
conditions (e.g., retarded, Feynman). The action is in- 
variant under infinitesimal coordinate transformations on 
the background spacetimc. 

We also need to expand the point particle action in 



powers of L. Using (2.6 1 we find 



S pp [z,g + h] 



f 1/2 

m I d\ {-g a ^z a z ) 



= 5(0) + 5(1) + 5(2) + . . . 
pp 1 pp 1 pp 1 

where the first two integration kernels are 



V 



a/3 



z<*zP 



pp (1) 



2m pi ( - g^z»z u ) 



1/2 



V 



z a zfz~>z s 



pp(2) 



4m 2 pl 



g^z^z" 



3/2 



(3.5) 

(3.6) 
(3.7) 



(n) 

and where S pp denotes the part of the action for the 
point particle containing terms proportional to n factors 

of huv. 



A. The closed-time-path effective action 

The construction of an effective field theory for the mo- 
tion of the effective point particle in a curved spacetimc 
begins with the CTP, or in-in, generating functional 



[ Vz» 


f VW V exp \ 


ICTP 


ICTP 1 



+ iS pp [z, g + h}- iS pp [z\ g + h'] 
+ i J dX(j^ - j'^z'n 

+ i [ d^g^iJ^-J^h'n], (3- 



where 

J CTP 

= j dzf J dz? J dz? £ Vz» £ Dz'" (.) (3.9) 

and likewise for the graviton CTP path integral. 

We find it more convenient to relabel the unprimcd and 
primed variables with a lowercase Latin index a, b, c, ... 
(from the beginning of the alphabet): 1 for an unprimed 
index and 2 for a primed index, respectively. We intro- 
duce the so-called CTP metric c a b that lowers and raises 
these indices where 



1 

C-ab \ q ! ( 



(3.10) 



For a current "contracted" with a scalar field, for exam- 
ple, the notation implies 

J Q $ a = c ab J a <t> b (3.11) 

= J!$! - J 2 * 2 = J 1 ® 1 - J 2 ® 2 (3.12) 

= J$-J'$'. (3.13) 

The Einstein-Hilbert action can be written as 

S[g + h a ] = S\g + hx]-S\g + h2] 

= S[g + h]-S[g + h']. (3.14) 

and similiarly for S pp . A comprehensive description of 
the CTP formalism using these notations can be found 
in [H|45]. In this compact notation the CTP generating 
functional takes on the form, 











[ Vz^ 


f Vh^ exp < 


iS[g + h a ] - 


F iSpp[z a: g 


ICTP . 


ICTP 







+ i / dXjffl + i / dW /2 «" 



(3.15) 



8 



Notice the similarity in appearance to the in-out gener- 
ating functional yet the substantive difference in physical 
consequences mentioned earlier. 

Calculating derivatives of the generating functional 
with respect to the external current J^„ generates time- 
ordered correlation functions of the quantum metric per- 
turbations hJi v 



s n z 



Z 5iJ^l Vl {xi) ■ ■ ■ 8iJ^ n 



(3.16) 



The CTP time-ordering operator T is defined so that un- 
primed operators are time-ordered, primed operators are 
anti-time-ordered and all primed operators are ordered 
to the left of unprimed operators. For example, the full 
graviton Feynman propagator is calculated from the gen- 
erating functional by 

iDf^ 5 (x,x') = (0,m\Thf(x)hJ S (x')\0,m) 

1 S 2 Z 
Z 5iJl p {x)5iJ] s {x>) 

(3.17) 



Notice that (3.16) describes the full graviton correla- 
tion functions and includes the effects from (nonlinear) 
particle-field interactions. Likewise, derivatives of the 
generating functional with respect to j!J generate time- 
ordered correlation functions of the worldline coordinates 
including the back-reaction from the metric perturba- 
tions. 

To guarantee a well-defined graviton propagator on the 
background spacetime we adopt the Faddeev-Popov [46] 
gauge-fixing procedure by introducing the action 



a 2 
s 9f = - m 



pi 



J d A xg l '*G a G a , 



(3.18) 



which is equivalent to picking the gauge G a [/i MIy ] ~ for 
the metric perturbations. (The « denotes weak equality 
in the sense of Dirac |47].) Since we will be dealing with 
tree-level interactions only there is no need to introduce 
ghost fields into the gravitational action. 

We choose the Lorenz gauge for the trace-reversed met- 
ric perturbations, defined as 



so that the gauge-fixing function is 







l a(3 



;0 



(3.19) 



(3.20) 



In this gauge, the kinetic term in (3.4 1 is 



d xg 



1/2 1 



1 



h- n h' a 



2h af> Rjfh 



7(5 



(3.21) 



which applies to both the h"^ and metric pertur- 
bations (equivalently, the unprimed and primed fields, 
respectively) . 

The generating functional can now be written as 



f Vz» fflq> J *S<°> [*«,]+* / 
J CTP { ' ' 



CTP 



oo oo ^ 

iS^[h a ] + iSg>[*a> K] • (3-22) 



n— 3 



n=l 



We have factored out the lowest order point particle con- 
tribution from the graviton path integral since it is inde- 
pendent of the metric perturbations. This is the limiting 
situation when the particle is regarded as a test object 
which produces no perturbations about the background 
metric (and thus its motion follows a geodesic of the back- 
ground geometry.) 

Perturbation theory in the in-in formalism is formu- 
lated similar to that in the in-out formalism. The fact 
that the metric perturbations couple linearly to the ex- 
ternal current J£ v implies that every occurrence of the 
field in (3.221 can be replaced by a functional derivative 



of the external current, 



K v {x) 



(3.23) 



This rule allows for the generating functional to be writ- 
ten as 



Z[f a , JH = I Vz» exp | iS<£ [z a ] + i I 
JcTP y J 

oo 

+ E 4? 

n=l 
oo 



0J af3 



n— 3 



0J ul3 



Zo[JH (3-24) 



since the interaction terms can now be taken outside of 
the graviton path integral. The quantity Zq is the free 
field generating functional for the metric perturbations 

z [JH =J c ^ p hflv ex P |^ (2) + 1 J d4 *5 1/2 wj 

(3.25) 

and is calculated by integrating the Gaussian giving 



Z a [JH = exp J - \ Jf ■ D^, s , • J^' (3.26) 
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where the free graviton two-point functions are 



D^^, s ,(x, x') 



-D 



-D + 



D S^' 8 ' ) (3-27) 



(momentarily leaving out the CTP tensor indices, also 
D c d — c ca CdbD ab ). Specifically, these are the Feyn- 
man (D F ) and Dyson propagators (D D ) and the posi- 
tive/negative frequency Wightman functions (D ± ). See 
Appendix [X] for their definitions, identities and useful 
relations. 

Upon defining the interaction Lagrangian as 



$ j a 



ft J a 



oo 

■£« 

n— 3 



S (r. 



0J uf3 



(3.28) 



we find that the generating functional can be written in 
the form 



, JH = I W exp I [z a ] + i I 

JCTP [ J 

<J 



+i I d A xC int 



• Z o[Juu\- 



(3.29) 



Notice that this is expressed as a certain functional 
derivative operator acting on a Gaussian functional of 
the external currents J£ v . 

Computing the Legendre transform of the generating 
functional with respect to the particle and graviton cur- 
rents gives the effective action 



r[(^),(C)l 



- / d 4 xgV 2 J"(kl- 



(3.30) 



The equations of motion for the expectation values of 
the worldlinc coordinates and the metric perturbations 
are then found by varying the effective action, 



= 



6T 



ST 

5(hsr> 



zi=z 2 ,h 1 =h 2 ,j a = J a =0 



Zl=Z2,hi=h 2 ,ja = Ja=Q 



(3.31) 
(3.32) 



In this paper, we are interested in deriving the equa- 
tions of motion for the particle including the radiation 
reaction force (self-force). We therefore focus our atten- 
tion on (3.31 ) and le ave th e field equations describing the 
gravitational waves ( |3.32[ ) for a future paper [T^j- This 
allows us to set J£ v = in the remainder of this paper. 



Before calculating the effective action we make a few 
remarks. A classical equation of motion does not ade- 
quately describe the particle's motion when the quantum 
mechanical fluctuations of the worldline are not negligi- 
ble. However, the particle's quantum trajectories can 
be decohered by interactions with the quantum fluctua- 
tions of the metric perturbations, or other matter fields 
present, resulting in a classical worldline for the particle. 
We call this the semiclassical limit (classical particle in 
a quantum field). The condition for the existence of a 
semiclassical limit and the appearance of stochastic be- 
havior are stated in the beginning of the Introduction 
section. See [TJ [2] and references therein for more de- 
tails of issues pertaining to a particle moving in an arbi- 
trary curved background. So as not to distract the reader 
with these issues it suffices to know that for a solar mass 
compact object (treated as a point particle) moving in 
the geometry of a 1O 5 M Schwarzschild black hole that 
the quantum mechanical fluctuations are of the order of 
10 _24 cm indicating that for EMM astrophysical objects 
the semi-classical limit is well-defined, as expected [2]. 



B. Power counting rules 

All of the terms following the kinetic term in 
(3.221 represent self interactions of the field and vari- 
ous particle-field interactions. Each of these interaction 
terms may be represented by a Feynman diagram. To 



write down all of the relevant diagrams that contribute 
to the effective action at a specific order in e we need to 
know how each of the interaction terms in (3.221 scale 



with e and 1Z. The scaling behavior that we develop here 
are called power counting rules and are essentially a gen- 
eralization of dimensional analysis. We first develop the 
power counting rules for the parameters of the effective 
field theory; we ignore for now the non-minimal point 
particle couplings in S pp (i.e., eg, cb, . . .). 

As discussed previously, the curvature scale 1Z de- 
scribes the length scale of temporal and spatial varia- 
tions of the curvature in the background geometry. This 
implies that each of the spacetime coordinates scale ac- 
cording to 



TZ. 



(3.33) 



From the kinetic term for the metric perturbations we 
deduce that if ~ 1 then 



I ^ R'lr (I) - R 2 lr 



and the metric perturbation scales with 1Z as 

1 



(3.34) 



(3.35) 



The particle-field interactions, indicated by the terms 

(n) 

Spp for n > 1, contain inverse powers of the Planck mass, 
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TABLE I: Power counting rules 



X ^ h[_LV 


L 


m 
m pl 


c( n ) 


S (n) 


* i 

1Z 










mlZ 


x/eL 




m p i. To power count these terms we form the ratio 

( m pi 



m p i 



(3.36) 



We remark that the factor mil is the scale of the (con- 
served) angular momentum for a test particle following 
a geodesic in the strong-field region of the background 
spacetime since 



L = mg a ^ a x^ ~ mil 



(3.37) 



where i 13 is the 4- velocity of the geodesic. We therefore 
find that 



m 
m p i 



eL. 



(3.38) 



(Neither does the quantity mil vanish if a test particle 
plunges radially toward the SMBH (for which the angu- 
lar momentum is zero) nor its value rely on the existence 
of any isometries of the background spacetime. In order 
to avoid problems of interpretation arising from scenarios 
where the particle's angular momentum vanishes one can 
identify mil as the leading order classical particle action, 
which is non-zero for all geodesic motions; see (3.39) be- 



low with n = [7S]. Regardless, we use the symbol L to 
denote the quantity mil in what follows.) 



The four scaling laws in (3.331, (3.351, (3.37) and (3.38) 



determine the power counting rules for identifying the 
appropriate Feynman diagrams that enter into the eval- 
uation of the effective action. See Table HI 

We now turn our attention to power counting the inter- 
action terms in (3.221. We consider first the diagram for 



the interaction of n gravitons with the effective particle 
worldline, as shown in Fig.([ljt). The curly line denotes 
a two-point function D^g , s , of the metric perturbation 
(i.e. of a graviton). The straight line denotes the ef- 
fective point particle. We remark that from the point 
of view of the gravitons, the particle acts as an external 
source that couples to the metric perturbations. As such, 
the straight line in Fig.([lji) does not represent the phys- 
ical propagation of the compact object, but acts as an 
external field, with the vertex denoting the particle-field 
interactions. 



(b) 



FIG. 1: Interaction vertices. Diagram (a) shows the interac- 
tion vertex for n gravitons, denoted by curly lines, coupling 
to a point particle, denoted by a straight line. Diagram (b) 
shows the self-interaction vertex of n gravitons. The labels 
ar, a,2, ■ ■ ■ and b are CTP indices, which take values of 1 and 
2. 



The power counting of n gravitons interacting with the 
effective particle is given by 



Fig. dl 



rn 



ft 



(3.39) 



The self-interaction vertices that result from the nonlin- 
earity of the Einstein-Hilbert action are given in Fig. (JTJd) . 
The power counting for the self-interaction of n gravitons 
gives 



Fig. Qd) = 



V 



1-9 



(3.40) 



From Table [I] we see that the power counting indicates 
that every type of interaction term involving any number 
of gravitons scales as L p where p < 1. 



C. Feynman rules and calculating the effective 
action 

We now turn to calculating the effective action from 
(3.301 for J^ v = 0. Standard quantum field theory argu- 
ments jlHl SSI |SU] demonstrate that the effective action 
is given by 



(sum of all 1PI 
connected diagrams 



(3.41) 



By "connected diagrams" we mean those contiguous di- 
agrams constructed using the Feynman rules for the in- 
teraction terms in (3.221. By "1PI connected diagrams" 
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we mean those connected diagrams that are one-particle- 
irreducible |49] 150] . However, we are only interested in 
those connected diagrams that contribute at the classi- 
cal level since the quantum corrections due to graviton 
loops on the motion of an astrophysical body are utterly 
negligible. Therefore, the effective point particle world- 
line is assumed to be totally decohered [75] and we will 
simply represent the expectation value of the worldline 
coordinates operators (z") by their semiclassical values 
z a 

■ 

A diagram with I graviton loops scales as L 1 ^ 1 , in units 
where h — 1. Therefore, classical processes correspond to 
those diagrams that scale linearly with L (i.e., tree- level 
diagrams) and provide the dominant contributions to the 
effective action so that 



iT[z a ] = iS$[z a ] + 



sum of all O(L) 
connected diagrams 



higher order 1PI 
graviton loop corrections 



(3.42) 



In this manner we construct a systematic method for 
computing the self-force equations order by order in e. 

The relationship between the connected diagrams, the 
interaction terms and the power counting is provided by 
the Feynman rules so that given a diagram at a given 
order in e we can translate these into mathematical ex- 
pressions. The Feynman rules are the following: 

1. A vertex represents the particle-field interac- 



tion iV 



pp («) 



or the graviton self-interaction 



iV^"' a2nplP2 as appropriate, 

2. Each endpoint and vertex is labeled by a CTP 
index, a, = {1,2}, such that a factor of 1 and 
(— l) Qi+1 is associated with each endpoint and ver- 
tex, respectively. 

3. Include a factor of the graviton two-point function 

, g ,(x,x') connecting vertices labeled by CTP 
indices a and b at spacetime points x and x' , 

4. Sum over all CTP indices, integrate over proper 
time for each particle-field vertex, and integrate 
over all spacetime for each graviton self-interaction 
vertex, 

5. Multiply by the appropriate symmetry factor, S. 

The symmetry factor S for a particular self- force diagram 
is found from the following expression 



1 



ni\ n 2 ! • • ■ n k \ 



(3.43) 



where k is the number of different types of vertices and n, 
is the multiplicity of a particular vertex appearing in the 
diagram such that n i equals the total number of 

vertices. We will show how these rules are implemented 
as we continue. 




FIG. 2: The diagram contributing to the first-order self-force 
described by the MST-QW equation. 



To derive the MST-QW self-force equation we only 
need those diagrams that contribute at 0(eL). From 
the Feynman rules for the interactions in Fig.© it fol- 
lows that there is only one such diagram at this order, 
which is shown in Fig.pl). Therefore, the effective action 
to first order in e is 



^^[z a ] = -im J dr a + Fig.© + 0(e 2 L) 



(3.44) 



where 



Fig.© 



1 \ / im 
2\) \2mpi 



2 

. i — i •* j 



a,b=l 



X z a Z a Dufa'S' ( Z a i z b ) z b Z b ■ (3.45) 



The symmetry factor here is 1/2! since there are rt\ = 2 

occurrences of the vertex Vpp in the diagram. 

Because the quantum fluctuations of the particle 
worldline are so strongly suppressed, it is reasonable to 
isolate the dominant contribution to the effective action 
and expand about zf — z% ■ Defining the difference and 
semi-sum (or averaged) coordinates 



z ?+ z 2 



(3.46) 
(3.47) 



and implementing the identities in Appendix [A] we find 
that such an expansion gives 



iT[z a ] = -im / dr z p g^ a a" + 



^4" fdr (dr'z^iT) 
*r& J J 



xw^[z«]V v D$ ys ,(zZ, z% )z\ 4 + 0(z 2 _) 

(3.48) 



where the 4-acceleration is 
<(r) = 



Dz% 

~dV' 



(3.49) 



r is the proper time associated with the worldline de- 
scribed by the semi-sum coordinates z" so that 



g a p{z+)z%z 



.0 



1 



(3.50) 



12 



and the tensor w' 1Q;/3l/ [z] is denned by 



+ vru . 



(3.51) 
(3.52) 



We remark that the CTP generating functional and 
the effective action provide causal dynamics for the ef- 
fective particle's motion since the retarded propa gator is 
the only two-point function that appears in (3.48 1. Inter- 



estingly, however, the difference in using the in-out versus 
the in-in formalisms is not overtly manifest at first order 
in e. We will see an explicit example of these different 
approaches when we calculate the second order self-force 
equations in the next paper. 



We observe that the retarded propagator in (3.481 is 
divergent when t' = t. In order to have a finite and 
well-behaved force on the compact object from the metric 
perturbations we will need to regularize this divergence 
and possibly renormalize the appropriate couplings of the 
theory. 



D. Regular ization of the leading order self- force 

The CS-EFT approach is founded in the theory of 
quantum fields in curved spacetimc [48, 51.. The renor- 
malization of divergences in this context has received 
much attention over the decades and a considerable body 
of techniques has been developed to remove these diver- 
gences in a systematic and self-consistent manner. We 
therefore find it natural to renormalize the divergence in 
(3.481 using these methods even if they are somewhat 



unfamiliar in classical gravitational problems. 

Of these approaches the method of dimensional regu- 
larization 52J is particularly useful. This regularization 
scheme preserves the general coordinate and gauge sym- 
metries of the theory but is also a natural choice to use 
within an effective field theory framework [9 , 53 54 . The 
reason for the latter can be seen from the problems that 
can develop when a simple cut-off regularization is used 
for the divergent integrals appearing, for example, in the 
Fermi four-point interaction theory of weak interactions. 
We refer the reader to [53] for the particular details of 
this theory. In this effective field theory the mass of the 
W-boson Myy is very large compared to the other masses 
(e.g., quarks) and typical momenta in the problem so that 
the action describing the low-energy theory is an expan- 
sion in powers of 1/Mw ■ As a result, using a momentum 
cut-off A one finds divergent diagrams at each order that 
scale like 



A 



w 



(3.53) 



where p is an integer. Since A represents the scale at 
which high energy physics becomes relevant it is natural 
to choose A ~ Mw- Therefore, all of the (power) diver- 
gences at each order contribute at 0(1) and the pertur- 
bative expansion in the EFT breaks down, unless when 



they are resummed in a particular manner 53J . This fea- 
ture does not occur with dimensional regularization since 
the dimensional parameter \x never shows up as a power 
H p but appears only in logarithms. This is true of any 
so-called mass- independent renormalization scheme |53J. 

The smearing prescription developed and implemented 
in [H H] to regularize the divergent part of the retarded 
propagator should not be used within our CS-EFT be- 
cause it is a mass-dependent regularization scheme. This 
is easily seen by looking at the shift in the mass of the 
electric point charge (See Eq.(3.17) and the discussion 
following in |2J in which the divergence is 



5m 



Ne 2 



, A 



(3.54) 



where N is an 0(1) number and A is a mass scale 
for smearing the divergence on the particle worldlinc. 
Matching onto the full (microscopic) theory implies that 
1 /A is of order the "classical radius" of the charge r ~ 
e 2 /m, which defines the length scale at which the vac- 
uum polarization induced by the charge's presence be- 
comes relevant. (Pair creation becomes important at this 
scale and marks roughly the length scale important for 
quantum elcctrodynamical processes [55].) Then 5m/ m 
is an O(l) contribution thereby causing the perturbative 
expansion to break down since the high-energy physics 
no longer provides a small correction to low-energy pro- 
cesses. Therefore, the smearing regularization of [TJ[5] is 
unsuitable for using within the CS-EFT framework. 

Instead, we use the dimensional regularization scheme 
below because of its practical case and because it allows 
for the effective field theory to be renormalized in a man- 
ner consistent with the associated perturbation series in 
e. Because we are applying a quantum field theoretical 
renormalization scheme to a classical gravitational prob- 
lem we provide below a somewhat pedagogical discussion 
of the steps in the regularization and renormalization of 
the effective action. 

The propagator is divergent in the limit when x' — > x 
and can be written as the sum of a regular and a divergent 
contribution 1481. 



D = D r 



D a 



(3.55) 



(we are temporarily dropping the spacetimc indices as 
well as the ret label). The renormalized propagator is 
defined by the finite remainder 



D' 



D-D 



div 



Pf{D) 



(3.56) 



where Pf stands for the pseudofunction of the quantity 
in parentheses and is well-behaved as a (regular) distri- 
bution when x' = x. Generically, quantum two-point 
functions and propagators are regarded as distributions 
and therefore only make sense when integrated against 
a test function. Let us therefore define j(X) to be such 
a testing function so that we can form the (divergent) 
integral 



1 = 



dX'D(z a ,z a )j(X') 



(3.57) 
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where we evaluate the propagator on the particle world- 
line z a (\). We refer the reader to Appendix [B] for our 
notations and definitions regarding distribution theory as 
well as to the excellent texts of [551 157] , 



The divergent integral in (3.571 can be written as 



lim /(e) + Hie) 



(3.58) 



where /(e) is the divergent part of the integral and H(e) 
is the finite part. These are related to the renormalized 
and divergent propagators through 



so that 



Fp dX'D(z a ,z a )j(X') 



«A — e poo \ 

lim I / + / dX' D(z a 1 z a ')j(X / ) 



(—0 



-oo J \+c 

\l T^BP 



dX'D^(z a ,z a )j(X'). 



(3.64) 



We may then write the worldlinc integral of the full prop- 
agator as 



lim 1(e) = I dX' D div (z a ,z a ')j(X') (3.59) 

' -oo 



H(0) = / dX' D ren (z a ,z a )j(X') (3.60) 

J — OO 

= Fp dX' D(z a ,z a ')j(X / ) (3.61) 



where Fp denotes the finite part of the divergent integral, 
in the sense of Hadamard [55] . in (3.57). These relations 



follow from the fact that the renormalized propagator in 
(3.561 is a pseudo- function and, by definition, generates 



the finite part of (3.571 



In general, as we discuss in Appendix [B] the divergent 
part can be written in terms of its power divergent terms 
as well as powers of logarithms [55] |57] 



JV 



0=1 



(3.62) 



as e — > 0. Here e is related to a cut-off momentum, 
e ~ 1/A, and is used for illustrative purposes. To identify 
the divergent part of the integral in (3.57) we need a 
representation for D dlv that allows for a clear separation 
of the divergent parts from the finite terms. 

We use a momentum space representation for the 
graviton propagator that is initially introduced by Bunch 
and Parker for a scalar field in an arbitrary curved space- 
time in [51]. Keeping only those terms in D d%v that are 
divergent amounts to expanding the propagator up to a 
specific order (n th ) in derivatives of the background met- 
ric when using Ricmann normal coordinates. Throwing 
away all higher derivative terms in the expansion, which 
are ultraviolet finite as can be verified by power counting 
the momentum integrals, results in the divergent struc- 
ture shown in (3.62). We may therefore write 



r\div r^BP 



(3.63) 



where the superscript BP stands for the Bunch-Parker 
momentum space representation. 



Returning to (3.581, the finite part of the integral is 



defined via the pseudofunction of the propagator in ( 3.56 ) 



Fp dX'D(z a ,z a )j(X') 



— OO 

oo 



+ / dX'D^(z a ,z a ')j(X'). (3.65) 



Using these expressions, we find that the first order self- 
force in (3.481 is given by 



im 2 f°° 

Zm pl J-oo 

\Fp (°° dr' Dl%, 5 ,{zl,zi)z 

I J — OO 



+ / dT>Df? api , 5 ,(zl,zi)zXzi\ 

J — OO J 

(3.66) 

where we parameterize the worldline by the particle's 
proper time. 

Ignoring terms that are higher order in m, the diver- 
gent contribution that arises from the second term is 

/oo 
dr'Df^, 6 ,(zl,zi)zXz s ;. 
-OO 

(3.67) 

The particle worldline is a geodesic of the background 
spacetime at leading order so that the velocity at r' is 
related to that at proper time r through 



zii.r')=g\(z + (r').,z + {r))z^r) 



(3.68) 



where g 1 ^ is the bi- vector of parallel transport, which 
parallel transports a vector at z + {t) to another point 
z+(t') along the unique geodesic connecting these points, 
namely, the leading order worldlinc of the effective par- 



ticle's motion. The divergent integral (3.671 can then be 
written as 

/oo 
dr'Df^, 5 ,{z%,zi) 
-OO 

xg^{zi,z%)g s ' (T {zi,z%). (3.69) 
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The integrand is now a rank-4 tensor at z+(t) and a 
scalar at z + (r'), which we write as 



D 



BP (n a' 
(n)al3 1 S\ Z +^ Z + 



(3.70) 

In [36] we derive the momentum space representa- 
tion of the Feynman propagator for metric perturbations, 
which is found to be 



D 



d d k 



2 k s k { 



Jk-y 



a(cd)b 



VacRbsdt + ^ad^bsci 
^ db^csdt ^cd^dsbi 




3 k 6 
4(d - 3) 



where we have ignored those terms that fall off as k~ 5 
in the integrand since these give finite contributions in 
d — 4 dimensions, y denotes the Riemann normal coordi- 
nates (RNC) of z a (r') with respect to the origin at z a {r) 
and C re t is the contour appropriate for the retarded prop- 
agator. The hatted indices represent tensor components 
evaluated in Riemann normal coordinates. 



Using (3.71) the divergent integral (3.671 in RNC is 

d d k 



dr , 

oo Jc r J^Y l 

8 d 2 - 6d+ 10 k s k l 



o ik -v , 



1 d - 3 w r 



2 d - 2 fc 2 



(d-2f 



R dSbiV 



m(n - a) -b 



(3.72) 



where we have used z b z c z d R a ~u = 0. We provide the 
computational details to regularize I mn in Appendix [c] 
using dimensional regularization and demonstrate that it 
vanishes. 

The remaining finite part of Fig.© is therefore 



Fig.© 



im 

"Hi 



dr ztw^ v [z + }W u 



Fp I dr' D^% l5 ,{z%,zi)ziz 5 ;. (3.73) 



Notice that we do not need to renormalize any parame- 
ters of the theory at this order since dimensional regu- 
larization sets the power divergence to zero. 

Having regularized the leading order contribution to 
the effective action in Fig.© we now compute the equa- 
tions of motion from (3.731). The variational principle 



ST[z a ] 
Sz^(t) 



= 



(3.74) 



yields 



ma /i (r) 



2"$ 



[z a ]Vu 



Fp / dr'D r a %, s ,(z a ,z a )P z" . (3.75) 

J — oo 

Using the definition of Hadamard's finite part of the in- 



tegral in (3.641 and the fact that the retarded propagator 
is zero for r' > r we see that the finite part is given by 

Fp dT'D r * 5 ,{z<*,z<*')zi'z s ' 



= lim 



f C dr'D r a %, s ,(z a ,z a ')p'z s '. (3.76) 

J — oo 



Inserting this into the equations of motion [77] gives the 
equation for the self-force on the compact object moving 
in a vacuum background spacetime 



2mt 



-w^ v \z a ] 



lim f dT'V u D r a e ^, s ,(z a ,z a ')P'z 5 ', 

e ~ J — oo 



(3.77) 

which was originally derived by Mino, Sasaki and Tanaka 
[57] and by Quinn and Wald [55]. 



IV. EFFACEMENT PRINCIPLE FOR EMRIS 

While it is intuitive to expect finite size corrections 
to be negligibly small at the linear order self-force one 
may be concerned with such corrections at higher orders. 
Specifically, at what order in e is the motion of the com- 
pact object affected by induced tidal deformations? In 
this section we answer this question using coordinate in- 
variant arguments and demonstrate that such finite size 
effects from a small Schwarzschild black hole or neutron 
star moving in a background curved spacetime unam- 
biguously enter the self-force at 0(e 5 ) and as deviations 
from a point particle motion at 0(e i ). We also discuss 
the tidal deformations of a white dwarf star, which are 
more sensitive to the curvature of the background geom- 
etry. 



A. Non-spinning black holes and neutron stars 

To begin we write down the effective point particle ac- 
tion that includes all possible self-interaction terms con- 
sistent with general coordinate invariance and worldlinc 
reparameterization invariance, 



Spp [z] 



-rn 



dr + CE / dr £„„£'"' 



z_=0 



+c B JdrB^B^ 



(4.1) 
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FIG. 3: Graviton scattering off the background of a static and 
spherically symmetric extended body (e.g., a Schwarzschild 
black hole, a non-spinning neutron star). 



where we have already used a field-redefinition to remove 
those terms involving the Ricci curvature tensor. The 
terms containing the square of the Ricmann curvature 
(and higher powers) represent the influence of the finite 
size of the body as it moves through space. This is seen 
by noting that the equations of motion no longer have 
vanishing acceleration so that the effective point particle 
does not move along a geodesic of the spacetime, 



ma' 



l (r) = 2c E £ a0 £ a ^ 



(4.2) 



Such deviation from geodesic motion is typical of tidally- 
distorted bodies. 

The coefficients ce,b are parameters that depend upon 
the structure of the extended body. We must therefore 
match the effective point particle theory onto the full 
theory in order to encode this "microscopic" or "high- 
energy" structure onto the long wavelength effective the- 
ory. The matching procedure involves calculating (coor- 
dinate invariant) observables in both the effective theory 
and the full theory [78]. By expanding the observable 
of the full theory in the long wavelength limit, where 
the effective theory is applicable, we can simply read off 
the values of ce.b as well as any other coefficients in 
(4.1). We outline the steps for a matching calculation 



and perform an order of magnitude estimation of Ce,b 
for a spherically symmetric compact object. 

Below, we power count the scattering cross-section for 
graviton Compton scattering shown in Fig. (J3j) , which rep- 
resents the scattering of gravitational waves in the space- 
time of the isolated compact object. 

We first compute the cross-section in the effective point 
particle theory. The scattering amplitude is computed by 



expanding the terms in (4.1l proportional to ce and cb 



to second order in the metric perturbations, 



S pp [z] 



■■■ + c E j dr (e^E^ + 2£$£^ v 

+ £$£ {x) » v + 2£$£^ v + ---) 
+c B f dr (b$B^ + 2B$B<®'»' 



(4.3) 



where the superscript denotes the number of metric per- 

(2) 

turbations appearing in that function so that B}1 v is pro- 
portional to kaphas, for example. From the power count- 



ing rules developed in Section IIIB| wc find that the scat- 
tering amplitude associated with Fig. ^ scales as 



%A & 



ce,b 



in 



pi 



1 

TV 



(4.4) 



where the 1/7?. 2 comes from the two spacetime deriva- 
tives in the Riemann tensor and the & denotes "and a 
term that scales as" . While the cross-section includes 
contributions from other terms in the effective particle 
action it will contain one term proportional to c 2 E bi 



Ml 2 



-E.B 



"pl 



1 

TV 



&• 



(4.5) 



where the pp subscript indicates that this is the cross- 
section computed in the effective point particle theory. 

We turn now to the scattering cross-section in the full 
theory. A cross-section represents an effective scattering 
area and the only scale present in the full theory of the 
isolated compact object is set by its size r m ~ m/m^. It 
follows that 



a full 



f(-) 



(4.6) 



where / is a dimensionless function. In the long wave- 
length limit where r m /TZ -C 1 the cross-section will con- 
tain a term proportional to TZ~ S , 



Vfull 



■■hrl 



(4.7) 



Since quantities computed in the effective theory ought 
to match those computed in the long wavelength limit of 
the full theory we conclude that 



l pl l m 



(4.8) 



V 



upon identifying the 1Z 8 terms in both a pp and <r/„ 



Using (4.8 1 we can estimate the order in e that the non- 



minimal terms appearing in the effective point particle 
action (4.1 ) will affect the motion of the compact object. 



The first diagram that the finite size terms will contribute 
is proportional to ce,b and is shown in Fig.([4^). This de- 
scribes the deviation from the pure point particle motion 
experienced by the compact object due to the inclusion of 
the non-minimal couplings to the background spacetime. 
This diagram scales with e as 



Fig.gi) 



Ce,B 



dr 



1 

TV 



e 4 L 



(4.9) 



and enters at fourth order. 

This diagram does not couple to metric perturbations; 
it persists in the absence of gravitational radiation. As a 
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CE,B 



(a) 




CE,B 



(6) 



FIG. 4: Lowest order contributions to (a) deviation from 
(minimal) point particle motion due to the tidal deformations 
of the compact object and (b) the self-force from the interac- 
tion of gravitational radiation with these deformations. 



result, while Fig.|4^,) will affect the motion of the particle 
it is not a correction to the self-force. To find the order 
at which the tidal deformations affect the self-force we 
power count the diagram in Fig. (Bp) to find that 



Fig.gJ)) ~ c E>B dT Qj) ^-\/iZ~e 5 L. (4.10) 

Finite size effects therefore enter the self-force at fifth 
order in e. 



B. Non-spinning white dwarf stars 



In (4.8 1 we assume that the size of the compact object 
is of the same order as its mass, r m ~ m/m^. While this 
holds true for black holes and neutron stars it does not 
for white dwarf (WD) stars. WDs are thousands of times 
larger than their Schwarzschild radius and subsequently 
experience stronger tidal effects than a BH or NS with 
the same mass. In fact, the tides may be so severe that 
the WD is tidally disrupted at some point along its or- 
bit about the SMBH. As a result, we expect that finite 
size effects may be (numerically) enhanced and alter the 
WD's motion, possibly at a lower order in e. 

To see how this can arise we define the ratio of the 
compact object's radius to its mass by f co so that 



= fcoGm = 



fco to 
32n mj 



and the bookkeeping parameter e becomes 



' rn 



mi,TZ 



(4.11) 



(4.12) 



with m/m p i 
hole then fi,h 



y/eL/fco. If the compact object is a black 
= 2. For a neutron star having a mass of 



1.4M© and a radius between 10 and 16 km it follows that 
f ns varies between 4.8 and 7.7, respectively. 

Using the new expressions for r m and e it is not difficult 
to show that 



ce.b 



f - 

J CO 



m 



pi 



m pl r m 



(4.13) 



and therefore the leading order finite size effects from 
induced tides scales as 



Fig. @i) ~ U 



s 4 L. 



(4.14) 



Since e cx f co it follows that Fig.fik) scales as f\ Q . We 
remark that this observation agrees with that in [BO] 
where the authors observe that the induced quadrupolc 
moment of a neutron star scales as r^ s . 

The parameter e generally depends on some distance 
scale, r, set by the orbit of the compact object. As 
a result, e takes on different values at different orbital 
scales. Similarly, in the post-Newtonian expansion the 
small parameter (the relative velocity of the compact ob- 
jects v) depends upon the scale of the orbital separation 
r through the virial theorem, v(r) ~ (Gm/r) 1 ' 2 . 

To determine how the leading order finite size effects 
from Fig.(EJi) depend on the orbital scale r we consider 
the SMBH background to be a Schwarzschild black hole. 
It then follows from ( |4.12| and fll.l) that 



n 



2 3 / 4 3 1 / 4 / c , 



M \M 



-3/2 



(4.15) 



As the orbital scale decreases we observe that e increases. 

For a WD, e can be significantly larger than for a BH 
or NS with the same mass because f w d is typically much 
larger than both foh and f ns . In fact, f wc i can be so large 
that Fig.(EJi) may be numerically enhanced and impor- 
tant at orders in e lower than the naive fourth order 
estimation given in the previous section. Such enhance- 
ment is consistent with our physical intuition that a WD 
experiences stronger tidal deformations than a BH or NS 
with the same mass. 

To demonstrate this enhancement consider a simple 
example. Assume that the radius of the WD equals the 
Schwarzschild radius of the non-spinning SMBH it orbits. 
Choosing the radii to be 6000 km each it follows that 
the masses of the WD and SMBH are m w 0.6M q [6T1 
and M = 4OOOM0, respectively. For these values, (4.111 
indicates that 



6750, 



(4.16) 



which is much larger than both fbh and f ns , as expected. 

For this scenario, how do the leading order finite size 
effects of ( 4.14 1 change with the orbital scale, r, of the 
WD? Fig.(5ji] shows a log-log plot of f w d£ 4 as a function 
of the radial distance r/M from the SMBH. The dark 
line is f w d£ 4 while the dashed gray lines are plots of e s for 
s = Q, . . . , 4. As the orbital scale decreases, the dark line 
crosses several dashed gray lines indicating that Fig.|4^) 
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r/M 



FIG. 5: The effects from the finite size of a white dwarf star 
can be enhanced as it orbits in closer to the SMBH. The 
white dwarf seems to undergo some form of tidal disruption 
by either tidal disintegration (triangle and circle) or Roche 
lobe overflow (square). In either case, tidal disruption may 
be numerically equivalent to a second order process. 



is numerically enhanced from the naive e estimate for 
BH's and NS's to order e s because of the stronger WD 
tidal deformations. It is the scaling of the non-minimal 
couplings ce.b in (4.131 with the fifth power of f w d that 



is responsible for this enhancement of the leading order 
finite size effects. 

However, this enhancement can not proceed indefi- 
nitely. If the tidal deformations become strong enough to 
transfer material from the WD to the SMBH (via Roche 
lobe oveflow) or to disintegrate the WD at its Roche limit 
then the WD is tidally disrupted and can no longer be 
described effectively as a point particle. This indicates 
that our particular construction of a point mass effec- 
tive field theory will no longer satisfactorily describe the 
binary. 

Using Newtonian estimates we find the Roche limit 
for the rigid (fluid) WD to be near r w 24M (r ~ 
47M). These scales are represented by a triangle (cir- 
cle) in Fig.([5|. The WD's Roche lobe begins to overflow 
and transfer mass that accretes onto the SMBH when 
r w 40 M and is denoted with a square in Fig.(|5]). 
In all three cases the WD is tidally disrupted, a process 
that occurs near the s = 2 line. 

All of this suggests that finite size effects from induced 
tidal deformations can be enhanced for a WD until the 
star undergoes tidal disruption at 0(e 2 ). If one is inter- 
ested in calculating the second order self-force on a WD 
these tidal effects may have to be taken into account at 
some point during the binary's evolution to accurately 
determine the gravitational waveforms. 

If the SMBH mass is increased to 1O 5 M we estimate 
that tidal disruption occurs much closer to the SMBH's 
horizon. Describing the WD as an effective point par- 
ticle is therefore valid over much of its orbital evolu- 



tion. Increasing the mass of the SMBH further we find 
that WD tidal disruption occurs inside the event hori- 
zon and is therefore ignorable with respect to observables 
and processes outside of the SMBH. Therefore, for the 
SMBH masses relevant for LISA's bandwidth the tidal 
disruption of WDs is likely to be negligible except per- 
haps near the plunge and merger phases for a < 1O 5 M 
Schwarzschild SMBH. 

If the WD is tidally disrupted before it plunges into the 
SMBH, one can conceivably construct a new EFT that is 
valid at scales much larger than the orbital scale of the 
binary in which the SMBH, the WD and the mass trans- 
ferring to the SMBH are treated collectively as an effec- 
tive point particle. Given the complexity of such a sys- 
tem an EFT description would be very useful. However, 
the matching procedure to determine the large number 
of relevant non-minimal couplings is likely to be difficult 
given that high order intrinsic multipolcs will need to be 
included to accurately describe this system. 

We recapitulate our results from this section with 
an explicit statement of the Effacement Principle for 
non-spinning compact objects. 

Effacement Principle for EMRIs. Tidally induced 
moments will affect the motion of a compact object at 
0(s 4 ). 

For a white dwarf orbiting a Schwarzschild SMBH this 
effect may be numerically enhanced until the star under- 
goes tidal disruption at 0(e 2 ), which may be relevant 
when the SMBH mass is less than about 1O 5 M for the 
particular example discussed here. 



V. CONCLUSION 

We develop an effective field theory approach for sys- 
tematically deriving the self-force on a compact object 
moving in an arbitrary curved spacetime without the 
slow motion or weak field restrictions. The EFT is a 
realization of the open quantum system paradigm in sys- 
tems with a large scale separation such that the system's 
induced fluctuations from the backreaction of the coarse- 
grained quantum field is utterly negligible |8J . An initial 
value formulation of quantum field theory is adopted here 
using the closed-time-path (CTP) formalism for the in-in 
generating functional, which guarantees real and causal 
equations of motion for the compact object. As an il- 
lustration of the procedures involved in our approach we 
showed how to derive the MST-QW equation describing 
the (first order) self-force on a compact object. The CTP 
formalism is needed for a calculation of the second order 
self- force as will be shown in [TTj . 

We describe the compact object as an effective point 
particle that is capable of accounting for tidally induced 
finite size effects. In calculating the effective action we 
encounter ultraviolet divergences stemming from a point 
particle interacting with arbitrarily high frequency modes 
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of a graviton field. Using Hadamard's partie finie to iso- 
late the non-local finite part from the quasi-local diver- 
gences we are able to implement dimensional rcgulariza- 
tion within a (quasi-local) momentum space representa- 
tion for the graviton propagator [36 . As such, all power 
divergences can be immediately set to zero implying that 
only logarithmic divergences are relevant for renormaliz- 
ing the parameters of the theory. At first order, the ef- 
fective action has a power divergence and may therefore 
be trivially regularized using dimensional regularization. 

In the spirit of an Effacement Principle we find that the 
finite size of the compact object first affects its motion 
at 0(e 4 ) for a non-spinning black hole and neutron star. 
For a white dwarf star we deduce that such effects may 
be enhanced until the white dwarf is tidally disrupted at 
0(e 2 ) in which case the effective point particle descrip- 
tion, and in particular the effective field theory developed 
here, breaks down. One may conceivably construct a new 
effective field theory by treating the supermassive black 
hole, the white dwarf and the accreting mass as an effec- 
tive point particle possessing many relevant non-minimal 
couplings to the background geometry describing the in- 
trinsic moments of this composite object. 

The leading order finite size corrections cause a de- 
viation from the motion of a minimally coupled point 
particle that is not caused by interactions with gravitons 
but is due to the torques that develop on the tidally de- 
formed compact object. On the other hand, the self-force 
is affected by the induced moments of the compact object 
at 0(e 5 ). 

In summary, the EFT approach has at least two ma- 
jor advantages over the existing approaches: It provides 
a systematic procedure for carrying out a perturbative 
treatment, and an economical way to treat the ultravio- 
let divergences. Our CS-EFT improves on the PN-EFT 
introduced in |10j in that it is valid for a general curved 
spacetime and not limited to slow motion or weak field 
conditions. These will prove to be of special benefit for 
higher order self-force calculations. We will apply these 
steps to calculate the self- force at second order in e [TT] . 
the gravitational radiation emitted by EMRIs [T2] and 
the self-force on spinning compact objects [13 . 
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APPENDIX A: DEFINITIONS AND RELATIONS 
FOR THE QUANTUM TWO-POINT FUNCTIONS 

In this Appendix we collect some definitions, identities 
and relations for the quantum two-point functions that 
are relevant for this work. 

The positive and negative frequency Wightman func- 
tions are defined as 

(*>*') = (VWVM) ( A1 ) 

G^ YS >(x,x') = (h YS '(x')h af} (x)), (A2) 

respectively. The angled brackets represent the quantum 
expectation value so that 



(6) = Tr 0(53<)d 



(A3) 



and /3(Ej) is the density matrix of the quantum field given 
on a hypersurface E{ at constant coordinate time x° = ti . 

The Feynman, Dyson, Hadamard and commutator 
(also known as the Pauli- Jordan function or the causal 
function) two-point functions are, respectively, 

G^ YS ,(x,x') = (Th a p(x)hy S '(.x')) (A4) 

G° 0YS ,(x,x') = {T*h al3 (x)h YS >(x')) (A5) 

G^ rs ,(x,x') = ({h af3 (x),h Y5 ,(x')}) (A6) 

G% 0YS ,(x,x') = ([h^(x),h Y s'(x% (A7) 

where T is the time-ordering operator and T* is the anti- 
time-ordering operator. The field commutator is inde- 
pendent of the particular state used to evaluate it. Given 



the Wightman functions in (All and (A2| we write the 
above two-point functions in the form (ignoring the ten- 
sor indices from here on) 

G F {x,x') = 0(t - t')G+(x,x') + 6{t' - t)G-{x,x') 

(A8) 

G D (x,x') = 6{t' - t)G+{x,x') + 9(t - t')G_(x,x') 

(A9) 

G H (x,x') = G + {x,x') + G-(x,x') (A10) 
G c (x,x') = G+{x,x') - G-(x,x'). (All) 

From these we define the retarded and advanced propa- 
gators by 

- iG ret (x, x) = 6(t - t')G c (x, x') (A12) 
+iG adv (x, x') = 6(t' - t)G c (x, x'). (A13) 

In terms of the other two-point functions, these propa- 
gators satisfy the following useful identities 

— iG ret {x,x') = G F {x,x')-G_{x,x') (A14) 

= G + (x,x')-G D (x,x') (A15) 

iG adv (x,x') = G D {x,x')-G-(x,x') (A16) 

= G+(x,x')-G F (x,x'). (A17) 
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The Feynman, Dyson and Hadamard functions are not 
all independent since 

G H (x,x') = G F (x,x') + G D (x,x') (A18) 
= G + (x,x')+G-(x,x'). (A19) 

Under the interchange of x and x' the Feynman, Dyson 
and Hadamard two-point functions are symmetric, the 
commutator is anti-symmetric and 

G+{x,x') = G_(x',x) (A20) 
G ret (x,x') = G adv (x\x) (A21) 

suggesting that the Wightman functions and re- 
tarded/advanced propagators are a kind of transpose of 
each other. 



APPENDIX B: DISTRIBUTIONS, 
PSEUDOFUNCTIONS AND HADAMARD'S 
FINITE PART 

In this Appendix we present the basic structure, con- 
cepts and definitions of distribution theory that are rele- 
vant for this work. The reader is referred to [56, 57\ for 
more information. 

Consider the set of functions that are infinitely 
smooth C°° and have compact support on any finite in- 
terval. These functions, called testing or test functions, 
form a set T>. A junctional f is a mapping that asso- 
ciates a complex number to every testing function in T>. 
A distribution is a linear and continuous functional on 
the space of test functions D and is frequently denoted 
by the symbols (/, 0) and /. 

For a locally integrable function f(t) we can associate 
a natural distribution through the convergent integral 



dtf(t)<j>(t) 



for some testing function € D. Notice that we are us- 
ing the same symbol to denote both the distribution and 
the function that generates the distribution. This is an 
example of a regular distribution. All distributions that 
are not regular are singular distributions and will be our 
main concern in the rest of this Appendix. An example 
of a singular distribution is the well known delta func- 
tional 5. As this is not generated by a locally integrable 
function S(t) (even as the limit of a sequence of locally 
integrable functions |57p it must be a singular distribu- 
tion. 

Often, a singular distribution gives rise to a singular 
integral, which can be written in terms of its divergent 
and finite parts. For the purposes of clarity and illus- 
tration it is best to consider a simple example. Let us 
compute the integral 



0® 
t 



<It 



(B2) 



for 0(t) a testing function in D and 9(t) the step, or 
Heaviside, function. This integral is obviously divergent 
since 1/t is not a locally integrable function at the ori- 
gin. Nevertheless, we may extract the finite part (in the 
sense of Hadamard [58 ) of the integral by isolating the 
divergences from the finite terms. 
To this end we write 



ct>{t) = m+-tm 



(B3) 



where ip(t) is a continuous function for all t. Putting this 
into |B2| and integrating gives 



t '■ 



= lim 



0(0) log b - 0(0) log e 



dtip(t) 
(B4) 



where we assume that the testing function 0(i) vanishes 
for t > b for some real number b. The finite part of (B2| 



is defined to be the remainder upon subtracting off the 
divergent contribution (s). In this case, dropping the log e 
term gives the finite part of the integral 



Fp 



dt 



tf>(t) 



t 



0(0) log b - 



dtijj(t) (B5) 



where the symbol Fp denotes the Hadamard finite part of 
the integral. Therefore, the divergent part of the integral 
is given by —0(0) log e. 

A distribution that generates the finite part of the in- 
tegral is called a pseudofunction, which we now calculate 
for this example. Inserting (B3) into the finite part (B5l 



gives 
Jo 



t 



lim 



3 °dt^ + 0(O)loge 



(Bl) Since 



0(0)= / dt6(t)cf>(t) = (5,<f>) 



. (B6) 



(B7) 



it follows that the finite part can be written as an integral 
of a distribution with a testing function 



Fp 



dt 



0(i) 



lim 



dt 



1 



■*(*)loge 



which defines the pseudo- function, 



Pf d ^ = 6 f + S(t) lim log, 

t t 6^0 + 



<Kt), (B8) 



(B9) 



Therefore, the finite part of the integral generates a 
pseudo-function (a regular distribution) that yields a fi- 
nite value when integrated with a testing function 

r dtPf e ^<t>{t)=F P r dt 9 ^^). (bio) 

J -co £ J —oo t 
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Quite generally, the value that the distribution assigns 
to a testing function will have a divergent part consisting 
of both power divergences and powers of logarthmically 
diverging terms so that 



N M 

p— 1 p— 1 



(Bll) 



for some appropriate integers N, M. This form for 1(e) is 
related to the so-called Hadamard's ansatz [58] and ap- 
pears often in regularizing divergent quantities involving 
two-point functions of a quantum field in curved space- 
time HOED- 



APPENDIX C: DIMENSIONAL 
REGULARIZATION OF EFFECTIVE ACTION 

In this Appendix we give the explicit calculation of 
the divergent part of the O(e) effective action. In Rie- 
mann normal coordinates, y a describes the coordinate of 
a point x' relative to the origin at x and is defined in 
terms of a tetrad e° at x through 



y & = -ei(x)*<*(x,x'). 



(CI) 



Here a(x,x') is Synge's world function, which numeri- 
cally equals half the squared geodesic interval between x 
and x', and a a = a'- a is proportional to the vector at x 
that is tangent to the unique geodesic connecting x and 
x' . See [63 for further details. Let us equate the points x 
and x' with z!t(r) and z^(t'), respectively, on the leading 
order (i.e., geodesic) particle worldlinc. 

The arbitrariness of the tetrad at x implies that we 
may choose 



e° (r) = -u +a (r). 



(C2) 



Hence, the 



if we maintain the condition e°e? = 5%. 
particle's 4- velocity is orthogonal to the components of 
the tetrad in the spatial directions, e l a u^_ = where 
i = l,,..,d—l. Using these relations and (CI I, the 



Riemann normal coordinate representation of the point 



(t') on the geodesic worldlinc is given by 



/ = e«(r)<(r)(T-r')=4?(T-r') 
implying that 

k-y= k & y & = -fco(r-r'). 



(C3) 



(C4) 



Integrating the divergent integral in (3.721 over s 



t' — t gives a delta function forcing to vanish so that 



l (r) 



i d — 3 



2 d- 2 

8i d 2 - 6d + 10 
T (d-2) 2 



1 



k l W 



(C5) 

where k is the d—1 dimensional spatial momentum with 
J k = J d d ~ 1 k/(2Tr) d ~ 1 and i,j = 1, ... ,d — 1 are indices 
for the spatial directions. 



The momentum integral in the first line of ( C5 1 can be 
integrated by giving a small mass m g to the graviton so 
that 



1 



(k 2 



2 7 r( d - 1 )/ 2 



Ud-2 



dk 



(2^)^ 1 r(^i)7 (k 2 + ml) 



2 ) , o, d-l 



1-d 



i r (a i 2 f 2 
( 47r )(d-l)/2 rxo) ^ TO s j 



(C6) 
(C7) 



for some positive integer a. Strictly speaking, the in- 
tegral in ( |C6[ ) does not necessarily converge for d = 4. 
However, by analytically continuing to other values for 
d we find that the integral can be made to converge. In 
this way, the divergence is rcnormalized via the analytic 
continuation and a finite result is obtained upon letting 
m g — > and then choosing d — 4. For a — 1 and d = 4— e 
the integral is 



1 



k 2 



47T 



0(e). 



(C8) 



In the second line of ( C5 1 we find that 

khi 



R 



1 _ B 

k 4 — ^60 



k 4 ' 



(C9) 



which vanishes because the background spacetime is vac- 
uous. 

Therefore, in the limit that the graviton mass goes to 
zero we find that ( C8 1 vanishes and 



I^(r) = e^I mn (r) = 



(CIO) 



as claimed. 
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